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SUMMARY 

This report presents a method for rapidly estimating the overall forces 
and moments at supercritical speeds, below drag divergence, of transport con- 
figurations with supercritical wings. The method is also used for estimating 
the rolling moments due to the deflection of wing trailing-edge controls. This 
analysis is based on a vortex -lattice technique modified to approximate the 
effects of wing thickness and boundary-layer-induced camber. Comparisons 
between the results of this method and experiment indicate reasonably good 
correlation of the lift, pitching moment, and rolling moment. The method 
required much less storage and run time to compute solutions over an angle- 
of-attack range than presently available transonic nonlinear methods require 
for a single angle-of-attack solution. 


INTRODUCTION 

The development of the NASA supercritical airfoil has led to aircraft con- 
figurations which have demonstrated significant increases in drag divergence 
Mach number. (See refs. 1 to 4 . ) As part of this development effort, the 
results of an experimental study to determine the effects of wing trailing-edge 
control-surface deflections on the static transonic aerodynamic characteristics 
of a transport configuration with a supercritical wing have been reported in 
reference 5. 

In recent years several three-dimensional transonic theories have been 
developed such as the methods due to Boppe (ref. 6) , Bailey and Ballhaus 
(refs. 7 and 8), and Jameson (refs. 9, 10, and 11). These methods require 
large amounts of computer run time and storage compared with methods based on 
linear theory. Linear theory, of course, will not account for the nonlinear 
characteristics of the flow; however, linear theory can frequently be used to 
estimate certain aerodynamic quantities in the transonic range. Therefore, the 
present study was undertaken to determine whether linear theory could be used 
as a rapid method of estimating overall forces and moments (as opposed to 
detailed pressure distributions) of transport configurations at supercritical 
speeds. The method includes approximate corrections for the effects of wing 
thickness and boundary-layer displacement thickness. The analysis has been 
performed on the transonic transport configurations reported in reference 5 
utilizing the vortex-lattice method of reference 12. In addition to calculat- 
ing lift and pitching-moment characteristics, calculations were made to deter- 
mine whether linear theory could be used to estimate the rolling moments due 
to wing trailing-edge control deflections. Because the experimental rolling- 
moment characteristics reported in reference 5 were obtained by deflecting con- 
trol surfaces on only one wing, the theoretical method of reference 12 was par- 
ticularly well suited for this study since calculations with this method may be 
performed on configurations which have an asymmetric geometry. 





The method of reference 12 typically requires 65000g storage and 17 sec- 
onds of execution time on the Control Data Corporation Cyber 1 75 system operat- 
ing under NOS 1.1. This calculation would determine the force and moment 
characteristics over an entire angle-of -attack range. For the same computer 
system, the method of reference 9 would require approximately 240000g storage 
and on the order of 90 minutes of computer time to calculate converged pressure 
distributions for only one angle of attack. 

SYMBOLS 

The International System of Units, with the U.S. Customary Units presented 
in parenthesis, is used for the physical quantities found in this paper. Cal- 
culations were made in U.S. Customary Units. The data presented in this report 
are referenced to the stability-axis system. The moment reference point was 
taken to be the quarter chord of the mean aerodynamic chord of the reference 
trapezoidal wing planform. 

b wing span, cm (in.) 

Cl lift coefficient, Lift/qJS 

C Lf0 lift coefficient at zero angle of attack 

ACj differential rolling-moment coefficient, defined to be the rolling- 

moment coefficient for a negative control deflection minus the 
rolling-moment coefficient for a positive control deflection of 
same magnitude for same control 

Cjj, pitching-moment coefficient. Pitching moment/q^Sc 

c m,o pitching-moment coefficient at zero lift 

^,1 pressure coefficient on wing lower surface, (pj - P^) /?„ 

c p,u pressure coefficient on wing upper surface, (p u - P^J/q^ 

^ c p = c p,l - c p,u 

c local chord, cm (in.) 

c mean aerodynamic chord, cm (in.) 

ith trailing-edge control surface (see fig. 1) 
free-stream Mach number 
p local static pressure. Pa (lbf/ft 2 ) 

p^ free-stream static pressure. Pa (lbf/ft 2 ) 
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free-stream dynamic pressure. Pa (lbf/ft^) 

S wing reference area of trapezoidal box extended to center line, 

m 2 (ft2) 

Uiocal local streamwise surface velocity that differs from due to 

thickness effects only, m/sec (ft/sec) 

free-stream velocity, m/sec (ft/sec) 

uj streamwise perturbation velocity due to thickness only at ith 

elemental panel, m/sec (ft/sec) 

w^ local downwash velocity at ith elemental panel, m/sec (ft/sec) 

x/c distance from wing leading edge divided by local chord 

y distance measured spanwise from plane of symmetry, cm (in.) 

a angle of attack, deg 

04 local angle of attack of ith elemental panel, rad 

Sf^ deflection of control surface F-^ (positive when trailing edge 

is down) , deg 

6* boundary-layer displacement thickness, cm (in.) 

n = 2y/b 

A sweep of wing quarter chord, deg 

Subscripts : 

1 lower surface of wing 

u upper surface of wing 

THEORETICAL METHOD 
Basic Approach 

Linear theory has been employed for the analysis of transport configura- 
tions with supercritical wings operating at transonic conditions. The analysis 
was accomplished by using the vortex-lattice method to model the configurations 
as planar lifting surfaces. The typical lattice used for this analysis is shown 
in figure 1 . 
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Because the experimental rolling moments were achieved by the deflection 
of trailing-edge control surfaces on only one wing, the geometry and flow field 
were asymmetric when these controls were deflected. Accordingly, the vortex- 
lattice method briefly described in reference 12 was chosen for the present 
study since this method can analyze asymmetric geometries and asymmetric flow 
fields. 

Preliminary analysis indicated that the inviscid thin-wing theory was suf- 
ficient for the prediction of both the slope of the lift curve and the slope of 
the pitching -moment curve over the angle-of-attack range of the data. However, 
the theory did not predict the actual values of the lift and pitching moment. 
Because of the linear nature of the lift and pitching-moment curves, it would, 
therefore, be sufficient to obtain a more exact solution at any angle of attack, 
such as the zero lift angle or, say, the angle for the design lift. For the 
present study, it was decided to correct the inviscid thin-wing theory at zero 
angle of attack. In order to obtain a more exact solution, the effects of wing 
thickness and boundary-layer-induced camber were included in the calculations. 
Details of these modifications are given in the next two sections. 


Thickness Modification 

The conventional method of solving the linear three-dimensional thickness 
problem employs a distribution of sources (ref. 13, for example). The elemental 
source strengths are determined by matrix inversion with the appropriate bound- 
ary conditions being applied. However, to avoid extensive modification to the 
vortex-lattice program chosen for this analysis and to be consistent with the 
order of accuracy of the present linear method, the technique presented in 
appendix A of reference 14 was chosen to account for the effects of wing thick- 
ness. This technique differs from conventional thickness methods because it 
accounts for the interaction of the thickness and the camber and results in a 
modification of the camber distribution, A general discussion of this inter- 
action may be found in reference 15. To aid in illustrating the application of 
this concept by the method of reference 14, the basic approach is reviewed. 

In reference 14, Rowe, et al., state that the local perturbation veloci- 
ties due to wing thickness can be accounted for in the linearized boundary 
conditions. In doing so, the following formulation is obtained: 

Wi t 

” = — ( a iUiocal) 

^OO ^OO 

where w^ represents the local downwash, represents the local angle of 

attack, and Ui oca i is defined as being the local steady streamwise velocity 
that differs from due to thickness effects only. For thin-wing theory, 

Uiocal i s identically equal to U^. However, accounting for finite thickness 
yields 

w i 

Uco 


= — “i^cc + U i> - a iP + — 
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where (uj/U^) represents the local streamwise perturbation velocity due to 
thickness only, nondimensionalized by the free-stream reference velocity U^. 
Values of the quantity 1 + (u^/U^) may be obtained from experimental results 
or by using the streamwise component of two-dimensional theoretical velocity 
distributions- In the case of two-dimensional theoretical solutions for the 
thickness effect, some improvements might be expected if simple sweep theory 
were used to compute the thickness-induced flow velocity normal to, say, the 
leading edge and the results transformed to the streamwise direction. However, 
reference 14 does not indicate the use of simple sweep theory and, since the 
intent of this modification was to implement and evaluate the simplified 
approach suggested in reference 14, simple sweep theory has not been used. 

The quantity 1 + (u^/U^) could be computed by the use of any two- 
dimensional solution method. For the present study, the conformal transfor- 
mation method of Theodorsen (refs. 16 and 17) was chosen for this calculation 
because solutions can be quickly obtained and the method was felt to be consis- 
tent with the order of accuracy of the basic approach. The calculations were 
made for the symmetrical thickness distribution of the supercritical airfoils 
modified by the Prandtl-Glauert transformation to account for subcritical com- 
pressibility effects. The thickness modification was implemented by first 
establishing the chordwise distributions of the modification near the tip and 
near the root of the wing and then linearly interpolating in the spanwise 
direction 1 along constant percent chord lines. No thickness corrections were 
applied for the fuselage. 


Boundary-Layer Modification 

The boundary-layer influence was accounted for on the wing by modifying 
the local angles of attack of the mean camber surface to include the induced 
camber effects of the boundary-layer displacement thickness. The additional 
thickness effects due to the boundary layer were not included in the analysis. 

No boundary-layer effects were computed for the fuselage. 

In order to estimate the effects of the boundary-layer displacement on the 
longitudinal aerodynamic characteristics, it was assumed that the boundary-layer 
displacement thickness could be calculated by use of a two-dimensional strip 
analysis. Initially, calculations were made by using the linear theory of ref- 
erence 19 which assumes that the flow is subsonic and that the boundary layer is 


^For the NACA four-digit family of wing sections, reference 18 indicates 
that values of the total surface velocity (expressed as a ratio to the free 
stream) for intermediate thickness ratios may be obtained approximately by 
linearly scaling the tabulated velocity ratios for the nearest thickness ratio. 
The chordwise distribution of velocity ratio predicted by the method of 
Theodorsen for the geometries of the present investigation was found to scale 
linearly with thickness ratio to the same order of accuracy that the velocity 
ratios for the NACA four-digit series did. The spanwise variation of the 
chordwise distribution of this velocity ratio (and, hence, of the modification 
to the boundary conditions) could, therefore, be approximately accounted for 
by linear spanwise interpolation and thus simplify the computational procedure. 
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incompressible. However, this theory did not provide the proper pressure dis- 
tribution for calculation of the boundary layer at transonic speeds, and the 
resulting boundary layers caused negligible changes in the lift and pitching 
moment at zero angle of attack. Hence, the two-dimensional transonic theory of 
Bauer, Garabedian, Korn, and Jameson (ref. 20) was used to compute the displace- 
ment thickness. 

In order to correct the lift and pitching moment at zero angle of attack, 
the boundary-layer calculations were made for a model angle of attack of 0°. 
Since converged solutions could not be obtained at a Mach number of 0.90, all 
calculations were made at a Mach number of 0.80. However, calculations indi- 
cated that the changes in displacement thickness over this Mach number range 
and for constant angle of attack (zero) are much less than the changes in dis- 
placement thickness over the angle-of-attack range investigated. Since the 
latter variation in displacement thickness is not being accounted for, it is 
not inconsistent to neglect the former variation. The experimental unit 
Reynolds number for all configurations was 9.84 ><10® per meter. All viscous 
calculations were based on this Reynolds number. 

Calculations based on simple sweep theory showed that the effect of sweep 
on the boundary-layer- induced camber was within the accuracy of the calcula- 
tions. Thus, to be consistent with the thickness calculation, the 2-D boundary- 
layer calculations were also made in the streamwise direction. 

Bauer, et al. (ref. 20) ignore the laminar boundary layer and assume that 
the turbulent boundary layer grows from zero thickness at a specified "transi- 
tion location." A correction for this assumption has been made in the present 
calculations by the estimation of the laminar Reynolds number at the model 
transition-strip location. The results for a flat-plate boundary layer with a 
zero pressure gradient in the streamwise direction were used. It was then 
assumed that the laminar and turbulent Reynolds numbers based on momentum 
thickness were equal at the model transition-strip location. Thus, a Reynolds 
number based on length could be computed for a turbulent boundary layer which 
began some distance ahead of the model transition-strip location. This Reynolds 
number was computed by the method of Sommer and Short (ref. 21) and provided a 
virtual origin of the turbulent boundary layer. The virtual origin was ahead 
of the model transition-strip location and was input as the "transition loca- 
tion" in the method of reference 20. 

The resultant boundary layers were found to vary spanwise approximately in 
a linear manner. For this reason, the same linear interpolation technique used 
for the thickness modification was used to compute the boundary-layer effects 
on the wing. Figure 2 shows the resultant upper and lower surface boundary- 
layer displacement thickness distributions used for interpolation on the 
A = 33° configuration. 
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COMPARISONS WITH EXPERIMENT 


Description of Models 

The configurations of this report had supercritical wing sections and area- 
ruled fuselages with vertical tails but no horizontal tails. The wings had 
quarter-chord sweep angles of 33°, 38.5°, and 42°. The planforms of these wings 
are shown in figure 3. The reference area for these wings was defined as the 
area created by extending the outboard leading- and trailing-edge line segments 
inboard to the wing center line; as a result, the reference planforms were trap- 
ezoidal. The moment reference point was taken to be the quarter chord of the 
mean aerodynamic chord of this reference trapezoidal wing planform. Only the 
33° wing had control surfaces. The controls were installed on the right wing 
only of the 33° configuration in order to simplify model construction and test- 
ing. It was anticipated that a reasonable estimate of the effects of differen- 
tial control deflection on both wings could be obtained by combining the effects 
of opposite deflections on the same wing. For this reason, the present compari- 
sons of rolling moment are based on a differential rolling-moment coefficient 
Ac i defined as the rolling-moment coefficient for a negative control deflection 
minus the rolling-moment coefficient for the corresponding positive control 
deflection. A detailed description of the 33° wing configuration as well as 
its experimental characteristics may be found in reference 5. The data for the 
other two configurations were obtained from unpublished results of wind-tunnel 
tests by the Langley Research Center. 


Longitudinal Aerodynamic Characteristics 

A comparison between the theoretical and experimental chordwise load dis- 
tributions at two span stations for the 33° wing configuration at a Mach number 
of 0.80 is presented in figure 4. In general, the theory underpredicts the 
experimental loads in the vicinity of the leading edge and overpredicts the 
experimental loads in the vicinity of the trailing edge. However, the dis- 
crepancy between theory and experiment is reduced by accounting for the effects 
of wing thickness and boundary-layer induced camber. 

A comparison between theoretical and experimental longitudinal aerodynamic 
characteristics of the 33° wing configuration with no control deflections at 
Mach numbers of 0.80 and 0.90 is presented in figure 5. The wing pressure dis- 
tribution data of reference 5 indicate that for these Mach numbers, shock waves 
and regions of supersonic flow are occurring on the wing. It can be seen that 
accounting for the thickness and viscous effects substantially improves the 
correlation between theory and experiment and results in a reasonably good pre- 
diction of the experimental characteristics. 

In order to assess the effect of wing sweep in the present method of com- 
puting wing lift and pitching moment, similar computations were made for the 


38.5° and 42° configurations. These results are shown in figures 6 and 7, 
respectively. It can be seen that similar improvements in the calculation of 
and Cjj, were obtained by accounting for wing thickness and boundary-layer- 
induced camber. The discrepancies between the theory and the experiment are 
most probably due to the flat-plate representation of the fuselage, the use of 
linear theory at transonic speeds, and the approximate manner in which viscosity 
and wing thickness have been accounted for. 

Similar results have been obtained by Hess (ref. 13) who estimated the 
effects of viscosity on a swept, tapered, and untwisted wing with a symmetrical 
airfoil. At a low free-stream Mach number, Hess demonstrated that inclusion of 
the boundary-layer displacement thickness in the theory improved the agreement 
between the theoretical and experimental chordwise pressure distributions and 
the spanwise lift distribution. 


Rolling-Moment Characteristics 

To properly compute rolling-moment effectiveness, both thickness and vis- 
cous effects should be included. However, because of the difficulty of perform- 
ing the boundary-layer analysis for an airfoil with a deflected trailing edge, 
no viscous effects on control effectiveness were computed. 

Insufficient data for comparison with theory were available for control 
F-| . (See fig. 1 for control-surface locations.) Comparisons of the theoretical 
and experimental rolling-moment coefficients produced by controls F2 and F3 
with the model angle of attack equal to 3.5° (approximately the design lift 
condition) are presented in figures 8 and 9, respectively. For the lower Mach 
numbers and low deflection angles, good agreement was achieved between the 
linear thin-wing theory and the experimental results. As would be expected, at 
the higher deflection angles or at the higher Mach numbers, agreement between 
the inviscid theory and experiment deteriorated. 

Figures 10 and 11 summarize these rolling-moment results for controls F2 
and F3 as a function of Mach number. These figures more clearly show the 
effects of Mach number for the various control deflection angles. From refer- 
ence 5, the drag divergence Mach number at the design C L of 0.50 for this 
configuration, with controls undeflected, is 0.916. It can be concluded from 
figures 10 and 11 that below the drag divergence Mach number and near the 
design lift condition, a reasonably good estimate of the differential rolling 
moment can be made by use of linear thin-wing theory without wing thickness or 
boundary-layer effects. Inclusion of the wing thickness effects consistently 
resulted in a negative increment in the differential rolling-moment coefficient 
for controls F2 and F3 over the range of deflection angles shown. 

Rolling moments produced by control F4 were substantially overpredicted 
by theory (figs. 12 and 13). This large discrepancy between the theory and the 
experiment may be attributed to aeroelastic effects or trailing-edge flow 
separation. 

Figure 14 presents a comparison between theory and experiment for rolling 
moments produced with multiple control deflections near the design lift coeffi- 
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cient. As before, rolling moments are reasonably well predicted below the drag 
divergence Mach number with the thin-wing solution. At high Mach numbers, 
linear thin-wing theory consistently overpredicted the control effectiveness. 


CONCLUDING REMARKS 

A study has been conducted in order to determine the extent to which 
linear theory can be used as a rapid method of estimating the overall forces 
and moments of transport configurations with supercritical wings operating at 
transonic speeds. 

Below the drag divergence Mach number, linear thin-wing theory gave a good 
prediction of both the slope of the lift curve and the slope of the pitching- 
moment curve against lift. However, the lift coefficient at zero angle of 
attack C Lf0 and the pitching- moment coefficients at zero lift Cjj, 0 were not 
well predicted by the thin-wing potential-flow theory. 

It was found that for Mach numbers below drag divergence, the prediction of 
Cl o and Cjn f0 was greatly improved by inclusion of wing thickness effects 
which account for the interaction of thickness and camber. With the additional 
inclusion of two-dimensional boundary-layer-induced-camber effects, as deter- 
mined from a nonlinear transonic method, the theory gave a reasonably good pre- 
diction of the lift coefficient C L and pitching-moment coefficient C^, over 
the linear part of the angle-of-attack range. 

For Mach numbers below drag divergence, lift coefficients near the design 
condition, and control deflection angles up to ±15°, thin-wing potential-flow 
theory gave a reasonably good estimate of the rolling moments due to deflection 
of inboard controls. Including the wing geometric thickness effects in the 
calculations generally resulted in a less accurate rolling-moment prediction. 

The effectiveness of the outboard controls was greatly overestimated. 

The linear theory method (including thickness and viscous effects) was 
substantially faster and required much less storage than presently available 
transonic three-dimensional methods. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 16, 1978 
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Figure 7.- Effect of geometric thickness and boundary-layer-induced camber 
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Figure 13.- Comparison of theoretical and experimental differential 
rolling-moment coefficients produced by two different deflections 
of control F4. cl = 3.50°; A = 33°. 
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